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o\ : 

^>«^ , A unifying theory of the denaturation transition of DNA, driven by temperature T or induced 

by an external mechanical torque F is presented. Our model couples the hydrogen-bond opening 
and the untwisting of the helicoidal molecular structure. We show that denaturation corresponds 
to a first-order phase transition from B-DNA to d-DNA phases and that the coexistence region is 
■ naturally parametrized by the degree of supercoiling o. The denaturation free energy, the temper- 

ature dependence of the twist angle, the phase diagram in the T, F plane and isotherms in the a, F 
plane are calculated and show a good agreement with experimental data. 
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D ■ Denaturation of the DNA, due to its essential relevance to transcription processes has been the object of intensive 
works in the last decades. Experiments on dilute DNA solutions have provided evidence for the existence of a thermally 
T ' driven melting transition corresponding to the sudden opening of base pairs at a critical temperature Tm [Q. Later, 
\ following the work of Smith et al. ||^, micromanipulation techniques have been developed to study single-molecule 
behaviour under stress conditions and how structural transitions of DNA can be mechanically induced. While most 
single-molecule experiments have focused on stretching properties so far, the response of a DNA molecule to an 
. external torsional stress has been studied very recently [^|^, sheding some new light on denaturation From a 
^ ' biological point of view, torsional stress is indeed not unusual in the living cell and may strongly influence DNA 
functioning H,^. 



C/2 



I 



o 
o 



> 



On 



■ O For a straight DNA molecule with fixed ends, the degree of supercoiling a — (Tw — Twq) /Twq measures the 



twist Tw {i.e. the number of times the two strands of the DNA double-helix are intertwined) with respect to its 
counterpart Twq for an unconstrained linear molecule. In Strick et al. experiment a A DNA molecule, in 10 mM 
PB, is attached at one end to a surface and pulled and rotated by a magnetic bead at the other end. At stretching 
forces of ~ 0.5 pN, sufficient to eliminate plectonems by keeping the molecule straight, a torque induced transition 
' to a partially denaturated DNA is observed. Beyond a critical supercoiling Uc — —0.015 and an associated critical 
torque Fc ~ — 0.05ey/rad the twisted molecule separates into a pure B-DNA phase with a = ctc and denaturated 
regions with a = —1. Extra turns applied to the molecule increase the relative fraction of d-DNA with respect to 
B-DNA. 

■ In this letter, we provide a unifying understanding of both thermally and mechanically induced denaturation 

■ transitions. We show that denaturation can be described in the framework of first-order phase transitions with control 
parameters being the temperature and the external torque. This is in close analogy to the liquid-gas transition, where 
control parameters are the temperature and the pressure. Our theory gives a natural explanation to the BDNA-dDNA 
phases coexistence observed in single molecule experiments ||^. We give quantitative estimates for the denaturation 
free-energy AG, the temperature dependence of the average twist angle A{9)/AT, the critical supercoiling CTc and 
torque Fc at room temperature in good agreement with the experimental data. Furthermore the dependence of the 

■ critical torque as a function of the temperature is predicted. 

^ ' Our model reproduces the Watson-Crick double helix (B-DNA) as schematized fig 1. For each base pair (n — 
O ■ l;-- -1-^)1 we consider a polar coordinate system in the plane perpendicular to the helical axis and introduce the 
^ ' radius r„ and the angle ipn of the base pair The sugar phosphate backbone is made of rigid rods, the distance 
between adjacent bases on the same strand being fixed to i = 6.9A. The distance hn between base planes n — 1 and 
, y n is expressed in terms of the radii r„_i,r„ and the twist angle 0„ = ipn — (fn-i as 
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Ci \ hnir„,r„-i,9n) ^ y - rl - rl_^ + 2r„r„_i cos0„ . (1) 

The potential energy associated to a configuration of the degrees of freedom (r„,(^„) is the sum of the following 
nearest neighbor interactions. 

First, hydrogen bonds inside a given pair n are taken into account through the short-range Morse potential 

Fm(r„) = D (e-'^('-"--R) - 1)^ with R = lOl . Fixing a = G.^A'^ 0,|ll|, the width of the well amounts to Sa^^ ~ 
0.5 A in agreement with the order of magnitude of the relative motion of the hydrogen bonded bases [|2j. A base 
pair with diameter r > rd = R + G/a may be considered as open. The potential depth D, typically of the order of 
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O.leV Ipl,^^ depends on the base pair type (Adenine-Thymine (AT) or Guanine-Citosine (GC)) as well as on the ionic 
strength. 

Secondly, the shear force that opposes sliding motion of one base over another in the B-DNA conformation is 
accounted for by the stacking potential ys(r„,r„_i) — i<; e~^(''"+''"-i~2^) (r„ — r„_i)^. Due to the decrease of 
molecular packing with base pair opening, the shear prefactor is exponentially attenuated and becomes negligible 
beyond a distance ~ 5b^^ = lOA, which coincides with the diameter of a base pair ||l0| , [lT| , |l3| , ^ . 

Thirdly, an elastic energy V6(r„, r„_i, 9n) = K[hn — H]'^ is introduced to describe the vibrations of the molecule in 
the B phase. The helicoidal structure arises from H < L. in the rest configuration r„ = i? at T = OK, VJ, is minimum 
and zero for the twist angle 6'„ — 27r/10. Choosing H = 3 A, we recover at room temperature T = 298K the thermal 
averages ~ 3.4A and (6'„) ~ 27r/10.4 [Q. The above definition of Vb holds as long as the argument of the square 
root in (nl) is positive, that is if r„, r„_i, 0„ are compatible with rigid rods having length L. By imposing Vf, = cxd for 
negative arguments, unphysical values of rn,rn-i,9n are excluded. As the behaviour of a single strand (r > r^) is 
uniquely governed by this rigid rod condition, the model does not only describe vibrations of helicoidal B-DNA but 
is also appropriate for the description of the denaturated phase. 

As will be discussed later, the elastic constant K = O.QlAeV/ is determined to give back the torsional modulus 
C of B-DNA estimated to C = 860 ± lOOl [|l5|,0 at T = 298K. The parameters of the Morse potential D and of 
the stacking interaction E we have set to fit the melting temperature — 350K of the homogeneous Poly(dGdT)- 
Poly(dAdC)-DNA at 20mM Na~^ see inset of fig 3. This melting temperature coincides with the expected 
denaturation temperature of a heterogeneous DNA with a sequence GC/AT ratio equal to unity at lOmM Na^ [0, 
as the A-DNA in the experimental conditions of Q. Among all possible pairs of parameters (D, E) that correctly fit 
T,„, we have selected the pair {D = O.lGeV, E = AeV/A^) giving the largest prediction for AG, see inset of fig 2, that 
is in closest agreement with thermodynamical estimates of the denaturation free-energy. 

When the molecule is fixed at one end and subject to a torque F on the other extremity, an external potential 
Vr{On) = —T9n has to be included. A torque F > overtwists the molecule, while F < undertwists it. 

The configurational partition function at inverse temperature (3 can be calculated using the transfer integral method : 



dipN {R.iPN\T''\R,0) (2) 

As in the experimental conditions, the radii of the first and last base pairs are fixed to ri — rpj = R. The angle of 
the fixed extremity of the molecule is set to ipi — with no restriction whereas the last one ipN is not constrained. 
The transfer operator entries read < r,(p\T\r',ip' >= X{r,r') exp{—f3{Vb{r,r',9) + Vr{9))}x{(^) with X{r,r') — 
Vrr' exp{—f3{Vm (r)/2 + Vm {r')/2 + Vs (r,r'))}. The VtV factor in X comes from the integration of the kinetic term; 
x{d) = lii0<6 = (p — (p'<TT and otherwise to prevent any clockwise twist of the chain. At fixed r, r', the angular 
part of the transfer matrix T is translationally invariant in the angle variables ip, ip' and can be diagonalized through a 
Fourier transform. Thus, for each Fourier mode k we are left with an effective transfer matrix on the radius variables 
Tk{r,r')=Xir,r')Ykir,r') with 
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The only mode contributing to Zr is fc = once (pN has been integrated out in (0). The eigenvalues and eigenvectors 
of To will be denoted by A, and ijjq (r) respectively with Aq > A]^ > . . .. In the N oo limit, the free-energy 

density /^'"^ does not depend on the boundary conditions on ri and rjv and is simply given by = — fc^TlnAo'"^. 

Note that the above result can be straightforwardly extended to the case of a molecule with a fixed twist number 
Tw = N£, e.g. for circular DNA. Indeed, the twist density i and the torque F are thermodynamical conjugated 
variables and the free-energy at fixed twist number £ is the Legendre transform of fr- 

We have resorted to a Gauss-Legendre quadrature for numerical integrations over the range rmin = 9.7A < r < 
fmax- The Morse potential Vm increases exponentially with decreasing r < R and may be considered as infinite for 
r < 9.7 A pT[ |. The extrapolation procedure to rmax oo depends on the torque value F and will be discussed below. 
Using Kellog's iterative method , the eigenvalues Ag'"'' and associated eigenvectors ip^'' (r) have been obtained for 

— . . (T) 

q = 0, 1,2. Like a quantum mechanical wave function, (r) gives the probability amplitude of a base pair to be 
of radius r. Two quantities of interest are: the percentage of opened base pairs P = J°° dr \'>plf\r)\'^, the averaged 
twist angle (9) = -dfr/dV. 

Results for a freely swiveling molecule at room temperature are as follows. V'o'"~°^ entirely confined in the Morse 
potential well and describes a closed molecule. Conversely the following eigenfunctions . . . correspond 

to an open molecule: they extend up to r^nax and vanish for r < r^. They are indeed orthogonal to another family of 
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excited states that are confined in the Morse potential with much lower eigenvalues. The shape of the open states are 
strongly reminiscent of purely diffusive eigenfunctions, ip^q'~^\r) ~ sin{q'!r{r — rd)/{rmax — T^d)) leading to a continuous 
spectrum in the limit rmax oo. 

This observation can be understood as follows. For r,r' > rd, the transfer operator To(r, r') is compared fig. 2 to 
the exact conditional probability p{r, r') that the endpoint of a backbone rod of length L is located at distance r' from 



the vertical reference axis knowing that its other extremity lies at distance r 17|. For fixed r, To and p both diverge 
in r' = r ± L and are essentially flat in between. The flatness of To derive from the expression of V;,: a rigid rod with 
extremities lying in r, r' may always be oriented with some angle 9* (— > at large distances) at zero energetic cost 
Vb(r, r',0*) = 0. As a conclusion, our model can reproduce the purely entropic denaturated phase. 

As shown fig 2, at a critical temperature = 350K, Aq'""'^^ crosses the second largest eigenvalue and penetrates, 

as in a first-order- like transition the continuous spectrum. For T > T^, the bound state disappears and in fig 2 

becomes the eigenmode with largest eigenvalue Jl8| . The percentage of opened base pairs P exhibits an abrupt jump 
from to 1 at T„i, reproducing the UV absorbance vs. temperature experimental curve for Poly(dGdT)-Poly(dAdC)- 

DNA The difference AG between the free energy ff~^^ of the open state (g = 1) and the free energy of 
the close state {q — 0) gives the denaturation free energy at temperature T, see fig. 2. At T = 298K, we obtain AG ~ 
Q.022eV in good agreement with the free energy of the denaturation bubble formation AG ~ 0.025ey estimated in AT 
rich regions The thermal fiuctuations in the B-DNA phase lead to an undertwisting A (6') /AT ~ —1.4 10~'*rad/K 
which closely agrees with experimental measures A{6)/AT ~ — 1.7 10^''rad/K |p^ . 

The presence of an overtwisting (respectively undertwisting) torque F > (resp. F < 0) strongly affects fg , 

leaving almost unchanged the single strand free-energy f^''- The denaturation transition takes place at T,„(F) 

see the phase diagram shown in the inset of fig 3. We expect a critical point at a high temperature and large positive 

torque such that i/'i is centered on R ||2^ . 

The supercoiling, induced by a torque at a given temperature smaller than T„i(F = 0) = 350K, is the relative 
change of twist with respect to the value at zero torque in the B-DNA state, a(r) — {(9)r — (6')r=o)/(^)r=o- In fig 3, 
we have plotted the isotherms in the cr, F plane. Horizontals lines are critical coexistence regions between the B-DNA 
phase, on the left of the diagram and the denaturated phase on the right (with a = —1). The left steep line is found 
to define a linear relation between F and a : T = Kg {{9)r — {9)r=o)- The slope Kg does not vary with temperature 
over the range 298 K< T <350 K and is related to the torsional modulus G of B-DNA through G = Kg{hn) /{kBT^ 
[15|. The value of K appearing in the elastic potential Vb and given above was tuned to ensure that G = 860 A 
|15,^. At room temperature, critical coexistence between B-DNA and d-DNA arises at torque Fc = — 0.035eF/rad 
and supercoiling ctc — —0.01. These theoretical results are in good agreement with the values Fc = — 0.05ey/rad, 
CTc — —0.015 obtained experimentally 

We plan to combine the present model with existing elasticity theories of DNA [||,|6| to understand the influence 
of an external stretching force on the structural transition studied in this paper. It would also be interesting to see 
how the above results are modified in presence of a heterogeneous sequence. 
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FIG. 2. With respect to Vb{r = 60, r', 61 = 0) (dashed-dotted line, plotted as a function of r' in units of ksT for T = 298K), 
— lnT()(r = 60, r') (full liuc) is flat and comparable to — In p(r — 60, r') (dashed line) up to an irrelevant additive constant. 
Similar curves are obtained for any r > ra. Inset: Free-energies f^^'^^ (lower curve) and f^^^'^^ (upper curve) of B- and 
d-DNA. AG equals Qm2eV at T = 298K and vanishes at Tm = 350K. 
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FIG. 3. Isotherms in the a,T plane. For T = 298K, the critical supercoiling and torque are ctc = —0.01, Fc = —0.035eV/rad 
respectively. Inset : phase diagram in the T, F plane. At T = 298K, Fc = —0.035eV/rad while in the absence of torque (dotted 
line), Tm = 350K. 
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